Momentum injection in tokamak plasmas and transitions to reduced transport 
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The effect of momentum injection on the temperature gradient in tokamak plasmas is studied. A 
plausible scenario for transitions to reduced transport regimes is proposed. The transition happens 
when there is sufficient momentum input so that the velocity shear can suppress or reduce the 
turbulence. However, it is possible to drive too much velocity shear and rekindle the turbulent 
transport. The optimal level of momentum injection is determined. The reduction in transport is 
maximized in the regions of low or zero magnetic shear. 

PACS numbers: 52.25.Fi, 52.30.-q, 52.55.Fa 



00 



: 

Oh- 



. 

O 1 

ryj ■ 

: 



CO 
> 
CO 
CO 

p 

ON 

o 
o 



X 



Introduction. In this Letter, we study the effect of ve- 
locity shear on turbulent transport in tokamaks to answer 
two questions: (a) what is the optimal momentum input 
that minimizes radial energy transport? and (b) under 
what conditions do abrupt transitions to reduced trans- 
port regimes occur? Experimentally, tokamak plasmas 
can develop regions of reduced transport where the tem- 
perature gradient is much higher than the typical value 
for the same energy input [H, leading to more stored en- 
ergy and better performance at less cost. In large toka- 
maks, these Internal Transport Barriers (ITBs) are found 
in regimes with low magnetic shear and with a net mo- 
mentum input by neutral beams |2|, |3| . In previous work, 
flow shear [343] and the Shafranov shift [|[ have been pro- 
posed as causes for the transition to reduced transport. 
Here we highlight the physical influence of the velocity 
shear and momentum input. Employing basic properties 
of the turbulent transport deduced from new numerical 
results 0, 0] we show that there is an optimal level of 
momentum input, and prove that abrupt transitions to 
reduced transport are possible because the steady state 
transport equations have several solutions, allowing for 
bifurcations. We also obtain the conditions under which 
transitions can occur. 

State of numerical evidence. Refs. 0, 0] studied nu- 
merically the effect of flow shear on the turbulent ion 
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FIG. 1: Schematic dependence of the turbulent energy flux 
Qt on the velocity shear je and the temperature gradient k. 



radial energy flux Q t with finite @ and zero 7] mag- 
netic shear. In Fig. [T] we sketch the dependence of Q t 
on the dimensionless parameters k — R/Lt and je = 
(B P /B T ){R 2 /vu)\diu/drl where L T = \d{\nTi) / dr^ 1 is 
the scale of variation of the ion temperature, uj is the ro- 
tation rate, vu is the ion thermal speed, r and R are the 
minor and major radius, and Bt and Bp are the toroidal 
and poloidal magnetic field. The energy flux is normal- 
ized by the gyroBohm value Q 9 b — {pi/R) PiVu, with 
Pi and pi the ion pressure and gyroradius. The curves 
in Fig. Q] are generated by a simple analytic model cho- 
sen to approximate the zero magnetic shear results of 
0. For every k, there is a minimum Q t , and for suffi- 
ciently small k, this minimum is zero. For large 7^, the 
parallel velocity gradient drives an instability that rekin- 
dles the turbulence 0, [l(| • The dependence of Qt on n 
and je is qualitatively similar for finite magnetic shear 
6], but quantitatively there is a considerable difference: 
with zero magnetic shear, the minima in Q t are much 
smaller for the same k, and the region of 7^ for which 
the turbulence is suppressed is wider. 

The turbulent flux of toroidal angular momentum n t 
was also calculated in 0, 0] ■ It is also normalized by the 
gyroBohm value H g B = (pi/ R) 2 Rpi- The dependence 
of n t on k and je has a remarkable property: defining 
the normalized turbulent diffusivities as \t — Qt/x- and 
v t = Ht/[(BT/ Bp)je], the turbulent Prandtl number 
Pr t — v t Ixt was found to be ap pro ximately independent 
of k and je and of order unity [16j |. 

Graphical analysis. We analyze a plasma heated by 
neutral beams. Consider a flux surface that contains the 
volume inside which the energy and momentum are de- 
posited. The ratio of the injected momentum and energy 
fluxes is Hb/Qb ~ Cvu/Vb, where V& is the beam velocity, 
C is a geometrical constant dependent on the angle of the 
beams and n?, and Qb are normalized by the gyroBohm 
values. Thus, Hb/Qb is a constant that only depends on 
the characteristics of the beam. In experiments, Hb/Qb 
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is usually of the order of 0.1 

To determine k and 7^, we need to solve the equa- 
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FIG. 2: Energy flux Q vs. temperature gradient re for a con- 
stant ratio Tl/Q oi momentum and energy input. 



tions Q = Qt + Qn = Qb and II = Tl t + Tl n = Tl b |17|, 
where Q n = X« K an d n„ = v n (BT / Bp)^/E are the col- 



lisional neoclassical energy and momentum fluxes 12 1. 
The dimensionless diffusivities \ n and v n are propor- 
tional to the ion-ion collision frequency and depend on 
the magnetic field geometry. Importantly, the neoclassi- 
cal Prandtl number Pr n — v n /Xn ~ 0-1 is smaller than 
the turbulent Prandtl number Pr t ~ 1 (l2| . 

To describe the solutions to Q — Qt, and II — lib, we 
plot the curves of constant Tl/Q in a (k,Q) graph, as 
given in Fig. [2] 18]. The beam characteristics determine 
the curve of constant Tl/Q. Then, given Q, the temper- 
ature gradient k is easy to read off the graph. 

To understand Fig. [2j it is convenient to consider the 
(7b, k) parameter space and search for points of inter- 
section of curves of constant Q and Tl/Q. From the gen- 
eral shape of the constant k curves in Fig. [1] we infer 
the contours of constant Q in the (te, k) plane, shown in 
Fig.E^a). The transport is purely neoclassical for k below 
the critical value k c (7b) (this corresponds to the points 
in Fig. [T] where Qt vanishes; it is not the linear stability 
boundary [19j). At k < k c , the constant Q curves are 
horizontal because Q n is independent of 7^. At k > k c , 
since neoclassical transport is usually much smaller than 
turbulent transport, the constant Q curves are approx- 
imately the constant Q t curves. We stress that k c (ie) 
is the curve of Qt = 0, but it is not a curve of constant 
Q = Qt+Qn- In Fig. [Ha), we have exaggerated the 
difference. 

The curves of constant H/Q are also shown in Fig.[3la). 
For k < k c , the transport is neoclassical, and for k 3> K c , 
turbulence dominates. Therefore 



K = {U/Q)- 1 Pr n (B T /Bp) lE iom< Kc , (1) 
K=(U/Q)- 1 Pr t (B T /Bp) lE for K > Kc . (2) 

In both regimes, the curves of constant 11/ Q are straight 
lines passing through the origin. Since Pr t > Pr ni these 
lines are steeper in the turbulent than in the neoclassi- 
cal regime. To transit from the former to the latter, the 
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FIG. 3: (a) Curves of constant Q (dashed lines) and constant 
11/ Q (solid lines). The thin red line is the critical temperature 
gradient re c below which there is no turbulence, (b) Sketch of 
the intersection between a curve of given Tl/Q and curves of 
constant Q with Qi > Q2 > Qs- The black dash-dot line is 
Eq. (0. (c), (d) Similar sketches for higher Il/Q. 



curves of constant II/Q must approximately follow the 
curve k c {ie) because neoclassical and turbulent trans- 
port are comparable in its vicinity. In Fig. [3ja), the 
transition from Eq. ([T]) to Eq. (J2j is shown in detail by 
exaggerating the difference between the constant H/Q 
curve and k c (7_e)- This piece of the curve of constant 
Tl/Q is crucial for bifurcations. 

Fig. [2] was produced using Fig. I^a). The intersections 
of constant Q and constant Tl/Q curves for k 3> k c cor- 
respond to the high Q section of the curves in Fig. [2j and 
the intersections for k < n c form the neoclassical straight 
line. The region in between, where for each value of Q 
and II we can find several values of k, is examined in the 
(7b, k) space below in the section on bifurcations. 

Optimal momentum injection. In Fig. [21 it is clear that 
to maximize n, we need to operate on the red dashed line 
K max(Q) that corresponds to the the maxima in n at con- 
stant Q in Fig. (3Ja). However, once there, any increase 
in k achieved by increasing the energy input Q is small 
because turbulent transport is very stiff. Therefore, the 
optimal operation is at the maximum critical tempera- 
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FIG. 4: (a) Asymptotic approximation to the curves of con- 
stant Q (dashed lines) and constant Tl/Q (solid lines) in the 
transition region between neoclassical and turbulent regimes, 
(b) Graphical representation of condition © for type B so- 
lutions to exist, (c), (d) The cases in which type B solutions 
no longer exist for small and large 11/ Q. 



ture gradient, K c , max = « c (7£, ma x), given in Fig. [5] and 
Fig. 03a) as a red star. As a result, the optimal temper- 
ature gradient is n c ,ma,x, the optimal momentum input 
is U/Q = (B T /B P )Pr n ('yE,ma,x/Kc,m a ,x), and the optimal 
energy flux is Q = XnK Ciinax . 

Conditions for bifurcations. Transitions can happen 
only when there are several values of n and je for given 
values of Q and IT. The curves of constant Q and con- 
stant U/Q can intersect in multiple points, as exempli- 
fied by the thicker line in Fig. [2] Fig. ISJb) is a sketch 
of the curves in the (je, «) plane that correspond to this 
case. For Q = Qi, turbulence dominates and there is 
only one solution, A\. If we decrease the energy input 
to Q2, there are three solutions A 2 , B 2 and C2, where 
C2 is neoclassical. A jump from A% to C 2 reduces the 
transport and increases the temperature gradient. If we 
continue decreasing Q to Q3, the constant Q and 11/ Q 
curves become tangent and there are two solutions. For 
Q < *?3j there is only one solution, which is neoclassical. 

Increasing 11/ Q above the value in Fig. G^b) gives the 
curves in Fig.^c). The solution C with the largest k is 



not purely neoclassical because the large momentum in- 
put causes a largeparallel velocity gradient, which drives 
turbulence [1, 0, [tj, E3] ■ At even larger 11/ Q, the situa- 
tion is as in Fig. EJd), where there is only one solution 
for each Q and bifurcations are not possible. It is easy to 
see how these cases are reflected in the large U/Q curves 
of Fig. [2j We now discuss the conditions for several so- 
lutions to exist. 

Compare Figs.^b) and[3Jd). The existence of several 
solutions is determined by the slope of the piece of the 
curve of constant U/Q that transits between the neo- 
classical and turbulent regimes. We study that region 
to prove that to have several solutions and hence transi- 
tions, U/Q and Q must be within a domain determined 
by the shape of the curve k c {ie)- Near this critical curve, 
Qt - Xt(7E)[« - k c(ie)}, where \t = K c (dxt/dK)\ K=Kc 
and An — k — n c <C k c . We also assume that Ptt re- 
mains approximately constant even for k ~ k c . Then, 
expanding in Xn/Xt ^ 1 and ordering Pr t ~ Pr n ~ 1, 
we find the curves of constant Q and U/Q to be 



Akq(7e) = g 



Q _ Xnl^c 



Xt 



AKn/Qdfi) 



X n n c n/ Q - Pr n {B T /B P ){ lE /K c ) 
Xt Pr t {B T /B P ){ lE /Kc) - U/Q 



(3) 
•(4) 



We plot these approximate expressions in Fig. Eta). The 
expression for Akjt/q(7.e) is only valid for the transition 
region between the turbulent and neoclassical regimes, 
i.e., for jei < Je < JE2, where jei and jE2 are the 
intersections between the curve k c (7.e) and the lines 
© and ©, i.e., Pr t (B T / B P )[ lE i/ k c {iei)] = U/Q and 
Pr n (B T / Bp)[-f E 2/ Kc(jE2)} = U/Q. These points of in- 
tersection are marked as red squares in Fig. Eltb), and as 
red dash-dot lines in Fig. HJa). 

The solutions to Q = Qt and U/Q = U^/Qb are 
given by Akq^e) = Ak u/q {^ e )- In Fig. gja), we re- 
cast Fig. [3Jb) in terms of Ak. To have several solutions 
we need type B solutions that we define as intersections 
where Ak'q^e^) < Ak^^q^e.b)- Here ' denotes dif- 
ferentiation with respect to 7^. This condition gives 



k'c(7e,b) > 



1 UB P 
P~r~tQB~T 



«c(7. 



E,B 



JE,B 



K{ 1e ,b). (5) 



In Fig. HJb), the dashed line is k' c , and the solid line is 
K(je)- Condition (|S|) is never satisfied near jei because 
Ak' u /q — >• — 00 there. 

Condition (J5J) defines an interval 
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n n 

< Q < Q 



(6) 
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outside of which there is only one solution for each Q. 
Inequality (0 is plotted for U/Q = U/Q\ m and U/Q — 
U/Q\m in Figs.|U[c) and@|d), respectively. At both ^e ,m 
and 7b. m, K ' c = K. In addition, je. m — 1E2, whereas at 
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FIG. 5: (a) Critical temperature gradient k c for zero (solid 
line) and finite (dashed line) magnetic shear i based on [?| 
and 0. (b) Region in the (IT/Q, Q) space where abrupt tran- 
sitions can happen. 



1e,Mi h" — K' ■ Then, and 7_b,m are given by 



K c (jE,m) — -pj, 



k"{1e,m) 



7b, m 



3 r t 7£,m ' 
Je,m k 'c(Je,m) 



B.M J 



(7) 

(8) 



Once 7_B,m and 7s,m are known, IT/Q| TO and II/QIm 
can be obtained from k'J^e) = K(p/E), which leads 
to IL/Q = Pr t {B T lBp)n' c ( lE )[ lE / Kc { lE )] 2 . The lower 
limit IL/Q\ m appears because K(^ E ) is bounded by its 
values at lE2 , (Pr 2 n /Pr t )(B T /Bp)(IL/Q)- 1 , and at lE1 , 
Pr t (B T /B P )(n/Q)^ 1 , tending to infinity for IL/Q -> 0. 
The upper limit IT/(5|m arises because increasing II/Q 
shifts the interval jei < 1e < 7B2 towards higher values 
of 7b and eventually n' c becomes negative. 
For every IL/Q, there is also an interval in Q, 



Qm<Q <Q 



M; 



(9) 



for which multiple solutions exist. We show the Q limits 
for IL/Q = 0.04 in Fig. [2] The lower limit Q m is Q 3 and 
the upper limit Q m is the meeting point with neoclassical 
transport. 

In Fig.[5fb), we give the domain in the (IL/Q, Q) plane 
where for each Q and IL/Q there are several solutions for 
k and je- We do so for zero 0] and finite @ magnetic 
shear, whose k c (ie) curves are given in Fig. EJa). The 
derivative k' c is clearly smaller for finite magnetic shear 
because the velocity shear is less efficient in quenching 
the turbulence. As a result, as we expect from ([5]), the 
region with several solutions is smaller. Thus, transitions 
are more probable at zero or small magnetic shear. 

Conclusions. In tokamaks, the turbulent transport of 
energy and momentum satisfies two properties: (i) given 
the energy flux Q, the temperature gradient k has a 
maximum possible value achieved at a finite flow shear 
7_e; and (ii) the turbulent Prandtl number Pr t is ap- 
proximately constant. Using these properties we have 
found the optimal level of momentum input. In addi- 
tion, employing the fact that the neoclassical Prandtl 
number Pr n is smaller than Pr t , we have shown that 



transitions to reduced transport can occur in an inter- 
val LI/Q\ m < IL/Q < LI/Q\ M - Below II/Q| m the flow 
shear is not sufficient to suppress the turbulence, and 
above IL/Q\m the large parallel velocity gradient drives 
strong turbulence. For each IL/Q, since bifurcations oc- 
cur only when neoclassical and turbulent transport are 
comparable and the effective Prandtl number is not con- 
stant, Q must be in the interval Q m < Q < Qm- Below 
Q m , the transport is neoclassical or close to neoclassi- 
cal, and above Qm-, it is mainly turbulent. The region in 
(IL/Q, Q) space where transitions occur is wider for small 
magnetic shear, which is consistent with both numerical 
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and experimental [2[ indications. 
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